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Fig. 6. S.D.(a'^Cpdb) vs. CO2 injected on column for simulated 16-bit data. The line indicates the theoretical limit of precision as a 
fianction of injection size, as calculated from Eq. (12). 



At natural abundance, a standard deviation of appt = 
1.0 is approximately equivalent to S.D.(5'^Cpdb) 
= 1 .0%o. We used Eq. (12) to predict the standard de- 
viation as a function of injection size at 16-bit resolu- 
tion, compensating for an open split ratio of 8.4: 1 . The 
integration window, W, was assumed to be constant 
at 10 s, E = 5000, and /?n = 3 x 10^ f2. A plot of the 
calculated limits compared to the observed precision 
at 16 bits is shown in Fig. 6. There is good agreement 
between theory and experiment. The calculated preci- 
sion is within a factor of five of the observed precision 
for all injection sizes. More striking, the calculated 
precision is a "lower limit", as nearly all the measured 
precisions lie above the theoretical prediction. The 
biggest discrepancies occur for large injection sizes, 
where the effect of quantization error is minimized, 
and other sources of error (e.g. contaminants) may 
dominate. 

Eq, (12) can also be used to demonstrate that quan- 
tization error should be negligible for signals acquired 
with 24-bit digitizers. Eq. (12) predicts that only 
0.6 pmol of CO2 to the IRMS should be necessary 
to achieve a precision of 0.5%o if quantization error 
is the only limiting factor. However, counting statis- 
tics dictate that a minimum number of ions must be 
formed to achieve a specified precision to overcome 
the shot-noise limit. Merritt and Hayes [11] give this 
equation as: 



(2 X 106)(1 -f R)^ 



EmN.R 



(13) 



where a^ is the shot noise limited standard deviation, 
R is the natural abundance isotope ratio, E is the ion- 
ization efficiency, m is the moles of CO2, and A^a is 
Avogadro's number. Substituting E = 5000 and R = 
0.011, we find that 6 pmol of analyte is required to 
achieve S.D. < 0.5%o. Therefore, when high precision 
24-bit boards are used, the effect of quantization error 
is superseded by shot noise. 

Our theoretical treatment gives insight into why 
curve-fitting is less sensitive to quantization error than 
summation. In the summation methods we have dis- 
cussed, imprecision in a single data point chosen as the 
background is multiplied throughout the background 
correction; for a peak width of 7^ data points, the to- 
tal quantization error scales as N. In contrast, the al- 
gorithms used in curve-fitting minimize the sum of 
squares between the fit curve and every data point. In 
curve-fitting, the quantization noise for each individ- 
ual point is averaged over the entire curve; for N data 
points, the total quantization error scales as 7/'/^. For 
a peak width of 10 s and a sampling rate of 10 Hz, 
this translates into a 10-fold reduction of quantization 
error. 

The theoretical treatment we describe is appropri- 
ate for understanding the effects of quantization er- 
ror on summation integration methods that choose 
single points on either side of the peak to define a 
background. It does not examine the limits of other 
classes of data reduction methods. An obvious im- 
provement to the summation integration method would 
be to average n points on either side of the peak, 
which would increase the effective number of bits of 
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Fig. 7. i'^Cpdb vs. CO2 injected on column for (a) individual summation and (b) curve-fitting algorithms. Data was collected at 16 bits 
and lOHzon the APP2003. 



the background measurement by «'/^, While this ap- 
proach could work well for isothermal runs with con- 
stant background, it is much less suitable for complex 
GC-CIRMS chromatograms, where it is not obvious 
which points should be averaged; that is, which points 
represent pure background and do not contain chemi- 
cal noise or the tail ends of peaks. Ricci et al, observed 
that the averaging method gives slightly higher back- 
ground values then other corrections [3]. They also re- 
ported that the dynamic background correction (which 
uses single points) yielded improved (5-values over the 
averaging method. Thus, a method may be insensitive 
to quantization error, but may still give worse results 
due to other variables, 

3.3. Improving precision on a 16-bit IRMS 

To test the effectiveness of curve-fitting on 
GC-CIRMS data acquired by low precision digitiz- 



ers, we ran multiple CO2 injections on an APP2003 
using 16-bit digitizers, and otherwise in similar 
fashion to the work on the FMAT252. A plot of 
^ ^Cpdb versus injection size is shown for summation 
(Fig. 7a) and curve-fitting (Fig. 7b). Fig. 8 shows a 
plot of S.D.(5'''Cpdb) as a function of injection size 
for both integration methods. The theoretical limit 
on the summation method, calculated from Eq. (12), 
is shown in the same figure as a dashed line. The 
observed precision for the summation method agrees 
well with theoretical predictions; most of the data 
points lie just above the lower limit curve. Almost 
15 nmol of CO2 on column are necessary to achieve 
a precision of <0.3% using the summation method. 
Usmg curve-fitting, only 0.76 nmol are necessary to 
reach that level of precision, a 20-fold improvement. 
This is greater than the two-fold advantage seen by 
curve-fitting the 16-bit data from the FMAT 252. 
One possible explanation is that the APP2003 data is 



8 

o 

St ■ 

I 4 

c *^ ^ 

M A 

u> 1 




♦ Surnmallon 
p Gurve-fitting 



Calculated Lirnit for 
Summation 



0,1 



1 10 

COj irtjected on coiurnn (nmpi) 



100 



Fig. 8. S.D.(5'^Cpdb) vs. CO2 injected on column for runs on APP2003 with 16-bit boards. The dashed line indicates the theoretical limit 
of quantization error on precision as a fiinction of injection size, as calculated from Eq. (12). 
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affected primarily by bit noise, while the FMAT252 
data has other sources of noise that cannot be elim- 
inated by curve-fitting. The FMAT252 accepts a 
three-fold higher flow rate than the APP2003, so there 
is likely more chemical noise in the FMAT252 signal. 
The relative immunity to quantization error with 
curve-fitting permits the IRMS to be run at lower 
inlet flow rates, which effectively mcreases quantiza- 
tion noise by decreasing the number of steps between 
background and peak with relatively little influence 
on chemical noise. The advantages of lower inlet 
flow rates are a longer lifetime for the filament, and 
reduced need for pumps and pumping capacity. These 
benefits, plus the reduced need for expensive ADC 
boards, should make high-precision GC-CIRMS 
more amenable to portable and low-cost applica- 
tions. In prmciple, statistical considerations define the 
lower limits of flow rates. However, counting statis- 
tics dictate that S.D. = 0.5%o requires 6pmol of CO2 
to the source for a typical continuous flow IRMS 
{E = 5000), and GC-CIRMS applications usually 
work well above this limit. Thus, modestly lower res- 
olution and inlet flow rates should not significantly 
affect performance, so long as appropriate integration 
techniques are used. 



4. Conclusions 

Data reduction using curve-fitting is more robust 
than the conventional summation method in the pres- 
ence of even modest levels of quantization error. Us- 
ing data obtained on high precision digitizers, the 
curve-fitting algorithm required several-fold less CO2 
to reach benchmarks of high precision (S.D, = 0.3, 
0.6, and 1.0%o) at any of the three simulated board 
depths (12, 14, or 16 bits). The poor performance of 
the summation algorithm was particularly noticeable 
at the 12-bit resolution, where S.D. <1.0%o could not 
be reached even at the maximum injection size al- 
lowed by the dynamic range of the Faraday cups. We 



have derived an expression that describes the influ- 
ence of quantization noise on isotope ratios calculated 
from raw IRMS data, and shown that it accurately 
predicts the lower limit of precision. Our theoretical 
treatment assumes that quantization error is uncorre- 
cted between the m/z = 44 and 45 signals, and is 
appropriate for any data reduction algorithm that uses 
single points on either side of the peak to describe the 
background. 

Curve-fitting substantially improved precision on 
GC-CIRMS data collected by an instrument with 
16-bit digitizers. The summation algorithm required 
15imiol of CO2 on-column to achieve a precision 
of S.D. = 0.3%o, while curve-fitting required only 
0.76 nmol. Thus, IRMS with 16-bit ADC boards 
achieved high precision for less than 1 nmol of C on 
column, a common benchmark for GC-CIRMS appli- 
cations, despite using 16-bit ADC boards. Lower inlet 
flow rates, enabling reduced pumping requirements, 
may be an important advantage in some applications. 
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